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We study an alloy system where short-ranged, thermally- 
driven diffusion competes with externally imposed, finite- 
ranged, athermal atomic exchanges, as is the case in alloys 
under irradiation. Using a Cahn-Hilliard-type approach, we 
show that when the range of these exchanges exceeds a criti- 
cal value, labyrinthine concentration patterns at a mesoscopic 
scale can be stabilized. Furthermore, these steady-state pat- 
terns appear only for a window of the frequency of forced 
exchanges. Our results suggest that ion beams may provide a 
novel route to stabilize and tune the size of nanoscale struc- 
tural features in materials. 
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The spontaneous formation of steady-state patterns 
have been extensively observed in many equilibrium and 
nonequilibrium systems [0. While for equilibrium sys- 
tems (e.g. ferrofluids, block copolymer melts, etc.) pat- 
terning is a result of the competition between repulsive 
and attractive interactions of different length scales, in 
nonequilibrium systems (e.g. reaction-diffusion systems, 
etc.), steady-state patterning is often the result of the 
competition between several dynamical mechanisms. A 
conceptual connection between the two classes of systems 
can sometimes be realized with the construction of Lya- 
punov functionals and effective Hamiltonians, by which 
steady-state pattern formation in dynamical systems is 
interpreted as resulting from the competition between 
different types of effective interactions. 

The kinetic Ising-type model with competing dynam- 
ics, and its continuum mean field counterpart, are in- 
struments by which we hope to understand a whole class 
of nonequilibrium driven systems P,p[ , ranging from fast 
ionic conductors to alloys under irradiation. The main 
ingredient of this model is the competition between two 
dynamics: one one hand, a thermally-driven mechanism 
trying to bring the system to thermodynamical equilib- 
rium; on the other hand, externally imposed particle ex- 
changes of a nature essentially athermal. The usual at- 
tempt has been to express the steady state of the system 
in terms of effective Hamiltonians and effective thermo- 
dynamic potentials. P. Garrido, J. Marro, and collabo- 
rators Q, were able to derive effective Hamiltonians for 
several types of ID Ising models with competing dynam- 
ics. Z. Racz and collaborators [^j, studied the relation be- 
tween the range of the externally imposed exchanges and 
the range of the effective interactions. In the context of 
alloys under irradiation, using a kinetic Ising-type model, 



Vaks and Kamyshenko p] derived a formal expression for 
the steady state probability distribution in terms of ef- 
fective interactions, while from a continuum perspective, 
Martin M studied the corresponding dynamical phase di- 
agram by an effective free energy. The possibility of pat- 
terning as a result of the competing dynamics has not 
been considered in these works. However, in the limit- 
ing case of arbitrary length external exchanges, pattern- 
ing has been recently observed in mean field and Monte 
Carlo simulations [g[. In this situation, the coarsening of 
segregated phases (magnetic domains in the Ising case) 
saturates, leading to a steady-state labyrinthine pattern- 
ing at a mesoscopic length scale. This microstructure is 
rationalized in terms of a competition between the at- 
tractive nearest neighbors interactions, and a repulsive 
electrostatic-like effective interaction. These patterns do 
not appear if the external exchanges are short range, e.g. 
when they occur between nearest neighbors |^,0. The 
behavior difference between these two limiting regimes 
raises the question of whether there exists a critical value 
for the range of external exchanges for patterning to oc- 
cur. The main objective of this Letter is to address this 
question, which besides having its own theoretical inter- 
est, is relevant to alloys under irradiation. Indeed, forced 
relocation of atoms in displacement cascades may extend 
beyond nearest-neighbor distances, especially in the case 
of dense cascades or for open crystal structures ||ll[ . 

To render the problem more concrete, let us consider 
a binary alloy with a positive heat of mixing, under irra- 
diation. Each time an external particle collides with the 
solid, a local atomic rearrangement is produced. These 
rearrangements have a ballistic component that mixes 
the atoms regardless of their chemical identity, trying to 
bring the system to a random solid solution. Due to its 
local nature, the ballistic mixing will relocate atoms in 
a region of characteristic radius R. The case i? — > oo, 
or arbitrary-length ballistic exchanges, has already been 
studied: The macroscopic governing equation is identical 
to that describing a binary alloy undergoing a chemical 
reaction A ¥^ B [|2|, and the one describing a block 
copolymer (BCP) melt [|3J. From the studies of these 
systems, the physics of this case is well understood. In 
terms of the frequency of forced exchanges F, it has been 
shown that while high values bring the system to a ran- 
dom solid solution, there is a critical value below which 
the homogeneous concentration profile becomes unsta- 
ble towards phase separation. As in spinodal decompo- 
sition, enriched regions form and coarsen. However, the 



characteristic length of the domains /, instead of growing 
indefinitely towards a macroscopic phase separation, sat- 
urates at a mesoscopic scale, loo- For F values close to the 
critical value, the steady-state concentration profile has a 
sinusoidal-wave appearance, with diffuse interfaces, what 
is referred to as the weak-segregation regime. For smaller 
F values, the concentration profile presents sharper in- 
terfaces, with a square-wave-like appearance, what is 
referred to as the strong-segregation regime. In each 
regime, the characteristic length has been shown to fol- 
low a power law with the exchange frequency Zoo ^ F""^, 
with an exponent (j> of 1/4 for the weak and 1/3 for the 
strong-segregation regime |l3| . 

In the case of ballistic exchanges of a finite range R, 
there should still be a critical value F(i?) below which the 
system phase separates. The question is whether phase 
coexistence takes place at a macroscopic or mesoscopic 
scale. In principle, we can predict that if certain F equi- 
librates a wavelength I for arbitrary length exchanges, 
finite range exchanges must also generate patterns when 
R ^ I. It is, however, difficult to determine a priori 
the exact conditions under which coarsening will saturate 
as a function of R and F. In this Letter we show that 
given a certain R, there is an interval [Fi(i?),F2(-R)] for 
the stabilization of patterns. Above F2 the homogeneous 
concentration profile is stable, and below Fi coarsening 
continues with the time, with the system separating into 
macroscopic phases. The extent of this interval for pat- 
terning decreases with R, reaching a zero value at a crit- 
ical, nonzero value Re, when Fi(i?c) = ^2{Rc) = Fc- For 
ballistic mixing with a radius smaller than Re, pattern- 
ing is not possible. We also show that, given a mixing 
distance R, there is an upper bound for the wavelengths 
attainable as F — > T^ {R). These conclusions are in agree- 
ment with recent Kinetic Monte Carlo simulations of bi- 
nary alloys under finite-range ballistic exchanges ||l4| , [l5| . 

We study the problem using a Cahn-Hilliard-type de- 
scription of one-dimensional fronts, simulating the walls 
of the labyrinthine patterns, and we construct a varia- 
tional formulation to investigate the solution. The equa- 
tion describing the temporal evolution is composed of 
two terms, one for thermal diffusion and another one for 
ballistic mixing 0: 

'at ~ ^ 



(1) 



dt dt 

Here we have chosen to represent the concentration field 
by a globally conserved order parameter ip{x.) so that the 
homogeneous concentration profile (solid solution) corre- 
sponds to "0 = 0. In the previous equation, the first 
term is simply given by MV^(^), where F is the global 
free energy, and we have assumed a constant mobility M. 
As for the second term, we actually need to perform a 
derivation. For that purpose, let us consider first bal- 
listic mixing occurring one dimensionally between planes 
along a crystallographic direction. The rate of change of 



concentration in the plane i due to interchange of atoms 
with the planes labeled j is given by Q 
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j- = -rj2wj{A - t.+j) = "F(^, - m, (2) 



dt 



where Wj is a normalized weight function describing 
the distribution of ballistic exchange distances, and the 
brackets denote the corresponding (discrete) weighted 
spatial average. The extension to the continuum is im- 
mediate, and we write the governing equation as: 
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Wfl(x) is now a continuous function peaked around the 
origin with a width proportional to R, and the average 
denoted by the brackets is defined as: 



Wi 



Wfl(x — x')7/;(x')(ix'. 



(4) 



In the limit i? ^- 0, the ballistic term reduces to a Lapla- 
cian term, expressing diffusion. In limit i? ^ 00, we 
recover the governing equation for the case of arbitrary- 
length ballistic exchanges. 

In analogy to what was done in the case of arbitrary- 
length exchanges, we seek to find a Lyapunov functional 
for this problem, which we shall refer to as the free energy 
functional of the system [ p6[. This idea actually traces 
back to the work of LeiblerpTJ] , and Ohta and Kawasaki 
p8[ , on block copolymer melts. This functional is given 
hy E = F + jG, and it is built so as to determine the 
kinetics: ^ — MV^(|^). G is a new term describing 
effective interactions related to the ballistic term, and to 
simplify the notation, we use 7 = F/Af for the rest of 
this paper. 

For F, we use a Ginzburg-Landau free energy 



F= [{-A-iP^ 



+ Bi^^ + C\W^\^) d^. 



(5) 



while G is expressed by a self-interaction of the form 

G=iyyV^(x)g(x-x')V(x')rfx(ix', (6) 

with g a kernel satisfying 

V2g(x - x') = -(5(x - x') - wr{^ - x')). (7) 

To proceed further, at this point we need to make a 
choice of the weight function wr. A Yukawa- type po- 
tential has been proposed by N. Goldenfeld [H. This 
form fits the observed distribution distances of ballistic 
exchanges for crystals under irradiation pi| , while allow- 
ing us to handle part of the minimization problem ana- 
lytically. In one dimension, wji{u) = i?/2exp(— |u|/i?). 



Having stated the problem, we start by performing a 
stability analysis of Eq. ^. For small perturbations of 
the form e"*+*'^^ around a constant profile ip = 0, it is 
straightforward to find the dispersion relation: 



M 



2Ak^ - 2Ck^ 






(8) 



A series of plots of this dispersion relationship for differ- 
ent values of 7 is shown in Fig. |[ As in the case of arbi- 
trary length ballistic exchanges, there is a critical value 
72 below which the homogeneous concentration profile 
becomes unstable. Below this value, there is a window of 
k values {ki,k2) for which the homogeneous solution is 
locally unstable, suggesting a wavelength selection. For 
smaller values of 7, the dispersion relation resembles the 
one of spinodal decomposition, suggesting macroscopic 
phase separation. 

To confirm these predictions, we use a variational ap- 
proach, based in minimizing the free energy functional E. 
Let us consider first the weak-segregation regime, where 
the choice of wr allows us to solve the problem analyt- 
ically, and then consider the strong-segregation regime, 
where we need to appeal to a numerical treatment. 

In the weak-segregation regime, we can perform the 
minimization of E by considering a sine family of para- 
metric functions, i^i^x) — a sin(fca;). We obtain the energy 
per unit length: 
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4(1 + Pi2) 



(9) 



Minimization of this energy is performed analytically, 
and concentration patterning, indicated by solutions with 
nonzero values of k and a are found for an interval in 
7. In reality, the value of 71 predicted by this paramet- 
ric functions is an underestimation. Before that value is 
reached, the energy per unit length for the macroscop- 
ically separated system becomes lower than the energy 
per unit length for the sine profile. The crossover point 
determines the actual value of 71 in the weak segrega- 
tion regime approximation, where patterning occurs in 
the interval given by: 



2%/3C + J6[(3 - V6)Ai?2 _ (76 - 2)C] 



71 - V(- 

72 = (A + C/R'f/2C. 



3i?2 



(10) 



This interval shrinks to zero at a critical value of Re = 
\JC I A^ corresponding to 7c = 2A^ jC . As a consequence 
of 71 being determined by a crossover of energies, a tran- 
sition towards macroscopic phase separation occurs at a 
finite value of k. As a result, there is a bound in the 
wavelength of the patterns for a given K. 



In the strong-segregation regime, we need to improve 
over the approximation of sine waves for a proper eval- 
uation of 71 away from the critical point. To this 
purpose, we propose to minimize the free energy func- 
tional using the tanh-sine family of parametric functions 
"0(2;) — Q tanh(?7i/fcsin(fc2;)). The parameter m serves to 
change the wave profile continuously from a sinusoidal 
type, to a tanh-like type with sharp interfaces, matching 
the concentration profile of an equilibrium interface. 

It is easy to minimize first with respect to the param- 
eter a. We obtain the expression, for a > 0, 



kE)=- 



(Aei - Ce-i + 7£4)^ 
4Be2 



(11) 



where the e, quantities denote the energy per unit length 
associated with each one of the terms in the free energy 
functional. Here we can see that the values of 7 and 
k that minimize the free energy are independent of the 
parameter B, which only relates to the amplitude a.. 

The energy per unit length (the e terms) cannot be 
obtained analytically for these functions, so we proceed 
with a numerical strategy. The free energy per unit 
length is computed by numerical integration, and the ac- 
tual minimization is performed by means of the subrou- 
tine MNFB from NETLIB [|o|, based on a secant Hessian 
approximation. Figure shows k versus 7 plots for a se- 
ries of values of R. The physical parameters A and C are 
set to unity. From this plot we obtain the dependency 
of 71 and the corresponding wave vector k\ (related to 
the maximum attainable wavelength) as a function of R. 
Figure is a double-log plot of these quantities for large 
values of i?, showing a power law dependency. A fit of 
the data with the power laws 71 = pR~^ and fci = qR~" 
yield the quantities p = 2.26 ± 0.03, 6 = 3.039 ± 0.006 
and q = 0.54 ± 0.02, a = 1.03 ± 0.01. Still, an almost 
perfect fit is obtained by the simple laws: 71 = 21 R? and 
k\ = l/(2i?), as shown in the figure. The latter rela- 
tionship has the physical interpretation that R is the pa- 
rameter determining the maximum wavelength. Further- 
more, the combination of these two relationships yield 
the power law for the R ^ 00 case: k ex 7^/^. For values 
of A and C not equal to one, the corresponding depen- 
dencies can be derived by dimensional analysis, giving: 
71 = 2\J AC jR? and k\ = 1/(2_R). Figure ^ summarizes 
the steady-state regimes in the 7-i? space. Patterning 
may occur in the present model when the range of the 
forced exchanges exceeds a critical value, with a maxi- 
mum wavelength proportional to that range. In the case 
of alloys under irradiation, typical R values range from 2 
to 10 A, suggesting that patterns up to 100 A could be 
stabilized. Experimental work to test these predictions 
is under progress. 
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FIG. 1. Growth factor a;(fc) for small perturbations of wave 
vector k. 
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FIG. 2. Wave vector fc as a function of 7 for several val- 
ues of R. Dash-dotted lines mark the limits of patterning as 
predicted by the sine profiles in the weak-segregation regime. 
The dashed line correspond to the large R power-law fit. In- 
serts (a) and (b) show the concentration profile for _R=10 in 
the weak and strong segregation regimes respectively. 
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FIG. 3. Double logarithmic plots of k\ and 71 versus i?, for 

large R, showing a power law dependency. 
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FIG. 4. Steady-state regimes as a function of R and 7. 

Asymptotics for 71 are indicated by dashed lines, for R ^ Re 
(sine profiles) and R ^ Re (power-law fits). Inserts show 
cuts of 3D Kinetic Monte Carlo simulations (to be presented 
elsewhere) of an FCC A50B50 alloy with ballistic exchanges 
of range R ~ 5/\/2 lattice parameters. 



